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Onsager symmetry implies that a Brownian motor, driven by a temperature gradient, will also 
perform a refrigerator function upon loading. We analytically calculate the corresponding heat flow 
for an exactly solvable microscopic model and compare it with molecular dynamics simulations. 
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Cooling techniques have a significant impact not only 
on our everyday life but also on the advances in science. 
We have come a long way from the primitive method of 
evaporative cooling, which our body conveniently uses 
when we perspire, over the evaporation by expansion of 
a cooling liquid in domestic refrigerators, to high-tech 
methods including laser cooling, magnetic cooling, ra- 
diative cooling or quantum cooling. Temperature is a 
measure of thermal fluctuations but the latter are not 
directly observable in macroscopic systems. Recent ad- 
vances in nanotechnology and molecular biology however 
allow to manipulate and even manufacture constructions 
on a molecular scale, where thermal fluctuations can no 
longer be ignored. To run such machines, we could copy 
the mode of operation from their macroscopic counter- 
parts. An alternative, and arguably more promising ap- 
proach, would be to utilize thermal fluctuations rather 
than fighting them. A well documented example is the 
Brownian motor P, Q , which generates power through 
the rectification of thermal fluctuations. In this letter, we 
present a novel method of microscopic cooling based on a 
Brownian motor in which, almost paradoxically, thermal 
fluctuations themselves can be harnessed to reduce the 
thermal jitter in one part of the system. 

Our Brownian motor Q (see Fig. 1) consists of two 
parts, a triangle and a flat paddle, which are rigidly 
linked and move as a single entity along horizontal tracks. 
Its motion is induced, following Newton's laws, by the 
random collisions with the particles of the gas in which 
it is embedded. Due to the random nature of these colli- 
sions, one expects a resulting rather erratic motion of the 
motor. The question of interest is whether, due to the 
asymmetry of the triangle, this random motion is charac- 
terized by a nonzero average speed in a given direction. 
When both motor units reside in a single compartment 
with a gas (or liquid) at equilibrium, we argue that such 
a sustained motion is impossible because it would vio- 
late the second law of thermodynamics. The construc- 
tion would be a perpetuum mobile of the second kind (or 
Maxwell demon), because the motion could be used, for 
example, to lift a weight, hence extract work out of the 
single heat bath. This impossibility can more deeply be 



understood from the fact that a system at equilibrium has 
a perfect time-reversal symmetry. Any movie of such a 
system, or of any subpart of the system (like monitoring 
the motion of our motor) is statistically indistinguish- 
able when played forward or backward in time. Hence 
any form of systematic translation is impossible. 

Having concluded that no sustained motion occurs at 
equilibrium, we modify the construction by inserting a 
perfectly insulating wall between both motor units while 
not hampering their joint motion as illustrated in Fig. lb. 
A noncquilibrium state is now easily realized by assuming 
that the particles in the separate lower and upper com- 
partments are at a different temperature, respectively T\ 
and Ti- To intuitively understand what happens in this 
case, we first consider the extreme case T± > T2 = 0, 
corresponding to fully stationary particles (no thermal 
motion) in reservoir 2 (upper reservoir with the triangle 
unit). When a particle hits the paddle in reservoir 1, it 
transfers horizontal momentum to the motor. It will do 
so in random quantities and equally so from the left as 
from the right. This momentum will be dissipated as the 
triangle hits the stationary particles in reservoir 2. How- 
ever, it is clear that the resulting slowdown will be much 
stronger when the motion of the motor is towards the left, 
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FIG. 1: An asymmetric object on horizontal tracks, (a) Per- 
forms Brownian motion with no average systematic drift when 
in thermal equilibrium, (b) Sustained average systematic mo- 
tion when subject to a thermal gradient. 
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with the flat face of the triangle transferring large mo- 
mentum to the stationary particles. When the motion is 
towards the right, the inclined sides of the triangle trans- 
fer less momentum allowing it to move farther, in the 
same way as a sharp arrow can penetrate deeper. Wc 
conclude that an average systematic motion is expected 
pointing to the right in the direction of the sharp angle. 

The above handwaving arguments are confirmed by 
an analytic calculation, which is furthermore exact in 
the limit of dilute gases Q • The average velocity of the 
motor is found to be : 



V = PxPiQ - sin 2 i 
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Here m and M are the mass of gas particles and motor 
respectively and M ^> m is assumed. The densities of 
gas particles in reservoirs 1 and 2, p\ and p2, are assumed 
to be small. The shape of the motor units is defined by 
the vertical cross section L and the half apex angle of the 
triangle 8 (see Fig. la). In agreement with the above 
discussion, we find for T\ > T 2 a positive average velocity, 
i.e., a velocity to the right. At equilibrium, T\ — T 2 , the 
velocity is zero (and also for 9q = tt/2, in agreement with 
the fact that the triangle reduces to a flat paddle without 
preferred direction). Note furthermore that equilibrium 
is a point of flow reversal: for T\ < T 2 , the motor has a 
negative velocity. Somewhat surprisingly, it then moves 
towards the left, i.e., in a direction opposite to the sharp 
angle! 

The transition from a Brownian motor to a Brownian 
refrigerator can be made by invoking a general stabil- 
ity principle that appears in several branches of physics 
under different names: Newton's action-reaction law in 
mechanics, Lenz's law of magneto-induction in electro- 
magnetism and, relevant to the problem under consider- 
ation here, the Le Chatclier-Braun principle in thermo- 
dynamics. The latter states that an action on a system at 
equilibrium induces processes that attenuate or counter- 
act the original perturbation. This principle of negative 
feedback is clearly an expression of stability of the orig- 
inal equilibrium state. Now, consider again the above 
construction, but at equilibrium with T\ = T 2 . We per- 
turb this state by applying an external force F < on 
the motor, moving it to the left. According to the above 
principle, we expect that the system reacts by a process 
that induces an opposing motion to the right. In the 
above discussion we have identified such a process: we 
need to increase the temperature in reservoir 1 with re- 
spect to reservoir 2. This entails the appearance of a heat 
flow from reservoir 2 to reservoir 1, cooling down reser- 
voir 2. Similarly, an opposite heat flow from reservoir 1 
to reservoir 2 will arise when applying a force F > 0. 

The above discussion may appear wishful and does not 
reveal neither the mechanism nor the amplitude of the 



process by which this miraculous heat flow appears. For- 
tunately, a simple but basic argument from the theory of 
linear irreversible thermodynamics £| confirms the pre- 
diction and allows us to evaluate the strength of the phe- 
nomenon. We note that the Brownian motor, described 
above, is the result of a cross-effect: a thermal gradient 
AT = T\ — T 2 produces an average motor velocity V. 
In irreversible thermodynamics, one calls J\ — V a flow 
and X 2 = I/T2 — 1/T% a thermodynamic force 0. Since 
the flow is zero at equilibrium, i.e., in the absence of the 
thermodynamic force, we expect to find, for a small tem- 
perature difference AT, T x = T + AT/2, T 2 = T- AT/2, 
a linear relation between flow J\ and force X 2 = AT/T 2 , 
namely: Ji — Li 2 X 2 . The value of the coefficient L\ 2 is 
found from Eq QJ: 
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Cross-processes are no exception. Well-documented 
cases are the Seebeck, Thomson and Peltier effects, where 
temperature differences induce electrical currents or vice- 
versa. The important observation in our discussion is 
that Onsager Q has identified a profound symmetry re- 
lation between any cross-process and its mirror process, 
deriving from the time reversal symmetry of underlying 
microscopic dynamics. More precisely, for any relation 
J\ = Li 2 X 2 , there is a mirror relation J 2 = L 2 \X\ with 
an identical proportionality coefficient L 2 \ — L\ 2 . In our 
case, J 2 is the flow associated to the temperature gradi- 
ent X 2 , i.e., it is a heat flow Qi^ 2 from reservoir 1 to 
reservoir 2. The particle flow is associated with a differ- 
ence in the chemical potential. In our setting, the latter 
is produced by the application of an external mechan- 
ical force F. The corresponding thermodynamic force 
is Xi = F/T. The relation J 2 = L 21 X U with Eq. Q, 
implies that the heat flow J 2 = Qi^ 2 is given by: 
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This heat flow is quite large, namely of the order of the 
power delivered by the external force (thermal speed of 
the motor, y/ksT/M, times F) multipled by y/m/M. 
Furthermore, in agreement with our foregoing discussion 
based on stability principles, the direction of heat trans- 
fer depends on the direction of the force F, in such a way 
that it activates an opposing Brownian motor. 

While heat is transferred from one reservoir to another, 
we have to realize that the motion induced by the force 
will have another, quite familiar effect, namely dissipa- 
tion by the frictional force 7V, where 7 is a frictional 
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coefficient. The dissipation will heat up both reservoirs 
and thus it may forestall the cooling. This "Joule" en- 
ergy stems from the power input P = FV. But with the 
velocity itself proportional to the applied force V = F/~f, 
the dissipation is proportional to F 2 . The quadratic law 
conforms with the fact that the dissipation is independent 
of the direction of the applied force. We conclude that 
the cooling effect is dominant for small forces. The ap- 
plication of such forces along with the use of either small 
motor molecules (M small) or larger substrate particles 
(m large) will induce a significant cooling effect in one of 
the compartments. 

All the above results, the explicit value of the Onsager 
coefficients and in particular Onsager symmetry, are con- 
firmed by a direct microscopic calculation. The start- 
ing point for this derivation is the following Boltzmann- 
Master equation for the probability distribution P(V,t) 
for the speed of the motor: 

dV'[W(V\V')P(V',t) - W(V'\V)P(V,t)] (4) 

where W(V|V') is a transition probability per unit time 
for the motor to change velocity from V' to V due to 
the collisions with the gas particles. This approach is 
exact in the limit of dilute gases. Analytic expressions of 
the transition probability for general convex objects are 
availabe for the dilute gases |3j. 

Expanding Eq. @ in Taylor series up to second order 
in y^m/M, we find the following explicit result (to linear 
order in F and AT): 



Ji — LuXi + L12X2 
3-2 — L21X1 + L22X2 
with Ji, J2, X\ and X2 as defined before, and 
T k Blll2 T 2 
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The friction coefficients in each of the reservoir are given 
by 



(5) 
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and the total friction by 7 = 71 + 72. L12 is as given in 
Eq. (J2J) and Onsager symmetry, L12 = L21, is confirmed. 

The above result also supplies us with the explicit ex- 
pression for Joule heating in each reservoir: 
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The ratio of the heat transferred from reservoir 2 to 
reservoir 1 (<$2-»i = — <3i— 2) over the Joule heat dissi- 
pated in reservoir 2 is thus found to be 
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If this ratio is larger than one, reservoir 2 will cool down. 
The following condition ensues in terms of the applied 
small (negative) force: 



> F > -2(1 - sin 9 ) k B T — Lp x . 



(10) 



A similar condition applies for the range of small positive 
forces, inducing cooling in reservoir 1. 

Besides the issue of Joule heating, we need also ad- 
dress the fact that the Brownian motor itself conducts 
heat 0- From Eqs. JSJ) and © we obtain the Fourier 
law, L22A2 = kAT with k = L22/T 2 . Assuming there is 
no other conductivity leak involved, one concludes from 
Eqs. @, © and ©, that total heat flow (J 2 = 0) will 
vanish at a relative temperature gradient AT/T of the 
order of the ratio of the drift speed F/j of the motor 
over the thermal speed i/feT/m of the gas particles. 
Note that this result in independent of the mass M of 
the motor due to the fact that both heat currents, orig- 
inating from thermal fluctuations, share the same 1/M 
dependence. 

To test the above predictions in a more realistic model, 
we have performed extensive two-dimensional hard disk 
molecular dynamics simulations (diameter a — 1). The 
system comprizes two separate rectangles (size=1200 x 
300), with periodic boundary conditions in both direc- 
tions, each containing 3600 hard disks (density p\ — P2 = 
0.01). The motor unit is as represented in Fig. 1 (L = 10 
and 6>o = 10°). Averages are taken over 30000 realiza- 
tions. In Fig. 2, we show the evolution of the tempera- 
tures in both reservoirs, starting with T\ = T2 = 1 (di- 
mensionless units k B = m = 1 and M = 2 are used). We 
apply a weak force F = —0.007 to the motor. The cooling 
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FIG. 2: The motor of mass M = 2m is driven by a constant 
force satisfying Eq. 11UL namely F = —0.007. The temper- 
ature of two reservoirs, Ti (dashed line) and T2 (solid line) 
are initially equal to 1. Clearly the temperature of reservoir 2 
(T2) is decreasing, while reservoir 1 is heating up. Note that 
the average temperature T av = (Ti + Tz)/2 increases due to 
the Joule heat dissipated in both reservoirs. 
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FIG. 3: The rate of energy change in reservoir 2 as a func- 
tion of the mass of the motor upon application of a force 
F — —0.007. The energy change encompasses both the heat 
transfer from reservoir 2 into reservoir 1 as well as the Joule 
heating in reservoir 2. Cooling is clearly observed. 



effect in reservoir 2 is clearly observed. To compare the 
theoretical results with the simulations, we plot in Fig. 3 
the net heat flow injected into reservoir 2, Qi^2 + Qji- 
The agreement between theory and simulation is, as ex- 
pected, very good for large values of M/m, but even for 
small values of M/m the order of magnitude is correctly 
predicted by the theory. 

Moving beyond the academic discussion, we need to 
address a number of issues of more experimental or tech- 
nological relevance. Under which conditions could the 
above discussed heat flow be observed or applied? The 
explicit result given in Eq. ||3J| was derived for an ideal 
gas. Molecular dynamics carried out at larger densi- 
ties however show that the order of magnitude of the 
heat flow remains the same. We can therefore use the 
above result to estimate the effect, even when we have 
in mind motors operating in an aqueous or lipid envi- 
ronment. Consider a macromolecule embedded in lipid 
bilayers, under physiological conditions, T » 300 K. For 




FIG. 4: A rotational version of the proposed refrigerator, (a) 
A chiral rotor is placed across an insulating wall, (b) A chiral 
biological molecule will channel heat across a membrane when 
a small external torque is applied. 



a typical value of its mass M — 10 5 m (m being typi- 
cally the mass of a water molecule) and for a force of 0.1 
pN, Eq. J3J) predicts (pi = p 2 ) a heat flow of the order of 
5 x 10~ 17 J/sec. On the other hand, using a typical value 
for the drag coefficient of a protein 7 ~ lOOpN • s/m, the 
Joule heat is found to be of the same order, and both 
heat effects are comparable. Hence, a negative net heat 
flow of this order of magnitude will appear for F slightly 
smaller or M slightly larger. Such a flow will cool down 
an aqueous reservoir of diameter 0.1 pui by one degree in 
about one minute time. 

We finally note that an alternative construction with 
rotational instead of translational motion of the motor 
would appear easier to realize. One can imagine a chiral 
molecule protruding on both sides of a membrane that is 
set into rotational motion by some polarizing force, see 
Fig. 4. As we have revealed above, such a construction 
could in principle function as a Brownian refrigerator, 
cooling down one side and heating up the other. Putting 
such constructions in parallel would amplify the heat cur- 
rent. Putting them in series would allow to channel heat 
out of the system following a predesigned pathway. Aside 
from the possible application in molecular biology, the 
aforementioned advances in nanotechnology should al- 
low to construct carefully designed devices of the type 
described above. As a final comment, we stress that the 
Le Chatelier-Braun principle implies the appearance of 
a Brownian refrigerator under very general conditions, 
whenever an opposing thermal Brownian motor can be 
activated in response to an applied external or thermo- 
dynamic force. 

We thank Chris Jarzynski for helpfull discussions. 
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